Exact dynamical state of the exclusive queueing process with deterministic hopping 
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The exclusive queueing process (EQP) has recently been introduced as a model for the dynamics 
of queues which takes into account the spatial structure of the queue. It can be interpreted as a 
totally asymmetric exclusion process of varying length. Here we investigate the case of deterministic 
bulk hopping p = 1 which turns out to be one of the rare cases where exact nontrivial results for the 
dynamical properties can be obtained. Using a time-dependent matrix product form we calculate 
several dynamical properties, e.g. the density profile of the system. 
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I. INTRODUCTION 

The one-dimensional asymmetric exclusion process, 
which can be regarded as the prototypical stochastic in- 
teracting particle system [1] , has been intensively studied 
in view of its non-equilibrium properties [3J, exact solv- 
ability [3J 0] and applicability to practical problems [5] . 
The state space for the exclusion process is the set of 
configurations of particles (in other words, the exclusion 
process has a "spatial structure"), and each particle can 
hop to its nearest neighbor sites only if the target site is 
empty ( "excluded- volume effect"). 

On the other hand, the queueing process is one of the 
basic stochastic processes in the field of operations re- 
search [BHH] . In addition to its practical relevance it often 
appears as effective model, e.g. in all kinds of jamming 
phenomena. Usually the spatial structure of the queue 
is neglected, i.e. the queues are regarded as "compact". 
However, often this assumption is not justified, e.g. in 
pedestrian queues. Therefore, recently a queueing pro- 
cess with excluded- volume effect (exclusive queueing pro- 
cess, EQP) has been proposed [5HTT]. On a semi-infinite 
lattice, particles enter the system at the left site next to 
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FIG. 1: Exclusive queueing process (EPQ) with deterministic 
bulk hopping (p = I). 
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the leftmost occupied site and leave the system at the 
rightmost site. In the bulk the particles move according 
to the rules of the totally asymmetric exclusion process 
(TASEP), see Fig.[l] 

The EPQ can be interpreted as a TASEP with varying 
system length which allows to analyze its stationary-state 
properties [9HTT] . In a more recent paper [12] , dynamical 
properties of the EQP were analyzed. Especially, for the 
deterministic bulk hopping case, dynamical behaviors of 
the average system length and the average number of par- 
ticles were investigated exactly. In this paper we derive 
more detailed results for the dynamical properties. 

The stationary state of the TASEP and some of its 
generalizations have been solved by means of the matrix 
product ansatz in the recent two decades [J]. The appli- 
cation of the matrix product ansatz to the calculation of 
non-stationary states is quite challenging and has been 
achieved only in a few cases so far, see e.g. [T3HT7] . In this 
paper we will introduce a matrix product dynamical state 
for the EQP, providing an explicit representation for the 
matrices. We utilize it for calculating typical quantities 
both in queueing theory and exclusion processes. 

Here we define the EQP as a discrete-time Markov pro- 
cess on a semi-infinite chain where sites are labeled by 
natural numbers from right to left (Fig. [IJ. A new par- 
ticle enters the chain with probability a only at the left 
site next to the leftmost occupied site (j — L). If there 
is no particle on the chain, a new particle enters at the 
(fixed) rightmost site (j — 1) with probability a. Each 
particle on the chain necessarily hops to its right near- 
est neighbor site if it is empty, i.e. we consider the limit 
of deterministic bulk hopping (p = 1). A particle on 
the rightmost site leaves the system with probability /3. 
These transitions occur simultaneously within one time 
step, i.e. we apply the fully parallel update scheme. Since 
we restrict our consideration to the case of deterministic 
bulk hopping p = 1 (the so-called rule 184 cellular au- 
tomaton), the stochasticity of the model is due to only 
the injection and extraction probabilities a and j3. 

In [T0Tll2| . the phase diagram of the EQP was de- 
rived. The parameter space is divided into two regions 
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Sec. IV we determine the density profile and the particle 
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FIG. 2: Phase diagram for the EQP with p = 1. 



(Fig. jSJ), the convergent phase a < and the diver- 
gent phase a > In the convergent phase, the system 
approaches a stationary state which can be written in a 
matrix product form. On the other hand, in the diver- 
gent phase, a stationary state does not exist, and the 
average length of the system (L t ) and the average num- 
ber of particles (N t ) increase asymptotically linearly in 
time t. On the "critical line" a — both (L t ) and 

(Nt) exhibit diffusive behavior, i.e. they increase being 
asymptotically proportional to \fi. 

In this paper, we investigate the dynamical (i.e. time- 
dependent) properties of the EQP in more detail. In the 
next section we write down the dynamical state (solution 
to the master equation) in a matrix product form. Using 
this form we investigate the waiting time, which is one 



of the basic quantities in queueing theory, in Sec. Ill In 



current profile. Concluding remarks are given in Sec. [V] 
In Appendix we review results on the usual (i.e. without 
excluded-volume effect) discrete-time queueing process. 



II. EXACT DYNAMICAL STATE 

For each site j we define the state variable Tj = 1 or 
corresponding to being occupied or unoccupied, respec- 
tively. For simplicity we impose the initial condition that 
there is no particle in the system, i.e. an empty chain. 
The state space is 

S = {U}U{ff<-ff l |<6N l ffj G {1,10}} 
= {0, 1, 10, 11, 110, 101, 1010, 111, . . . } (1) 



which is a subset (S C S) of 

s = {0, 1} u {ir L _! • • • n \L - l e N,tj e {l, o}} . (2) 

The element corresponds to the state where there is no 
particle in the system. Note that, for p = 1, the sequence 
00 never appears if the system starts from the empty 
chain. For simplification, we do not write the infinite 
number of 0's left to the leftmost particle. We denote 
the probability of finding a state r G S at time t by 
Pt(r), and the initial condition is written as Po(0) = 1 
and P (t) = (r € S \ {0}). We denote the system 
"length" for the state r £ 5 or r e S by |r|, which is 
nothing but the position of the leftmost particle (Fig. [I]) . 
In particular we define |0| =0. 

For the generic choice of the parameters < a < 1 and 
< p < 1, the process is irreducible and non-periodic on 
S. The master equation is simply written as 



P i+1 (0) = (1 - a)/3P t (l) + (1 - a)P t (0), (3) 

P t+1 (l) = (1 - a)P t (10) + (1 - a)(l - P)P t (l) + aP t (0), (4) 

Pt+i(ul0) = (1 - a)pP t (lul) + (1 - a)£P t (10ul) + a/3P t (ul), (5) 

P t+ i(ul01) = (1 - a)P t (lul0) + (1 - a)Pt(lOtilO) + aP t (uW), (6) 

P t+X (ull) = (1 - a)(l - P)P t (lul) + (1 - a)(l - P)P t (Wul) + a(l - /3)P t (ul) (7) 



r 



for u € S. In particular, for u — we set 010 = 10, 
101 = 11, and so on. This simple form is due to the 
deterministic hopping p = 1. 

We derive an exact dynamical state, beginning with 
the factorization ansatz 

P t (r) = Q t (|r|)y(r). (8) 

The first part Q depends only on time and the system 
length, and the second part Y is independent of time and 



satisfies the following relations: 

r(«il01«a) =f3Y( Ul lu 2 ), (9) 

Y( Ul nu 2 )=(l- P)Y( Ul lu 2 ), (10) 

r( Ul io) =pY( Ul i), (ii) 

Y(l) =1, (12) 

F(0) =1. (13) 
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One can easily see that the solution to these relations is 
Y(t l •• - n) = - p)*#V\Ti=i}-L- n ( 14 ) 

for tl ■ ■ -Ti e 5\{0}. The relations l^-Q also have the 
following matrix product representation, which is more 
convenient later: 



Y(r L ---n) = (W\X TL ---X T1 \V) 



with 
X X = D 



1-/3 
VP o 



Xq — E — 
(W\=(ly/0), \V) 



o VP 



1 

VP 



(15) 

(16) 
(17) 



These are essentially the same matrices and vectors as 
for the matrix product stationary state for the EQP with 
p = 1 [H] . The first part Qt(L) gives the probability that 
the system length is i at time t since 

£ F(r)= J2 Y(t) = {W\D(D + E) l - 1 \V) = 1. 



|t|=£ 



r\=L 



(18) 



Note that we can replace S by 5 in the above equation 
thanks to E 2 = 0. Inserting the relations (|9|-(12) into 
the master equation ([3|-([7|, we obtain 

Q t+ i(0) =(1 - a)Q t (0) + 0(1 - a)Qt(l), (19) 
Qt+i(L) =aQ t (L - 1) + (1 - a)(l - /3)Q t (Z) 

+ (1 - a)j8Q t (i + 1). (20) 

These equations actually agree with Eqs. (65) and (66) in 
|12j which were derived in a different way. The solution 
to this recurrence formula with the initial condition 



?o(0) 



1. 



)(£) = 0(ie 



is given by [12] 



Qt(L) = C zt 



1 - A 

1-z 



-A 1 



with 



A = 



l-(l-a)(l-P)z-r 



2(1 -a)Pz 

r = v/[l - (1 - a) (I - /3)z] 2 - 4(1 - a)apz 2 



(21) 
(22) 

(23) 
(24) 



where C z tF(z) denotes the coefficient of z* in the Lau- 
rent series for the function F(z), that is C z tF(z) — 
§ 2iriz t + 1 F( z ) with a small anti-clockwise path enclos- 
ing the origin of the complex plane. The average system 
length (L t ) at time t is derived as [TZ] 



(£ i ) = ^LQ t (L) = C^- 



A 



L>0 



(1-*)(1-A) 

(«< lf?)> 
(a 



1+0 < 



(25) 



(26) 



^(a-p + ap)t (a>j^), 
for f — > oo. 

Inserting Eqs. (151 and (22 1 into Eqn. (j8j), we obtain 
the matrix product dynamical state 



Pi 1 



C z t - 



i't(Ti---Ti) = C7, 1 



1 - A 

1 - A 

l-z 



A i (W / |X T 



(27) 

•X Tl |F). (28) 



When a < (convergent phase), the matrix product 
dynamical state converges to the matrix product station- 
ary state pj] 



lim P t \ 



lim(l - A) 



[3 — a — a/3 



(29) 



t^oo " ' z ^r (3(1 - a) ' 

lim P t (T L ---n)= lim(l-A)A L (W\X TL ---X Tl \V) 

t— fOO Z— >1 



/3 — a — a/3 
P{l-a) 



(l-a)/3. 



(W\X TL ---X T1 \V). (30) 



III. WAITING TIME 

The waiting time is one of the most important quanti- 
ties in queueing theory, which corresponds to the number 
of time steps that a particle needs to leave the system af- 
ter entering the system. 

Before we derive the waiting time distribution we de- 
termine the distribution of the number N of particles in 
the system. In standard queueing theory N is always 
identical to the length L of the system since the queue 
has no internal structure. In the EQP we only know that, 
by definition, N can not be larger than L. The probabil- 
ity that the number of particles is N = at time t is, of 
course, equal to Qt(0). For N € N, we find 
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P A+B {N) = P A (N) + P B (N) = J2 Pt{T L ---T l ) = C zt l Y -^{W\(kD + k 2 DE) N \V) 

t l --t 1 £S: 

#{j|T i= l}=W 

A(l-A)(l + W K 
[1-z) 

where P t A (N) [ resp. P B (N)] is the probability of finding N particles in the system and the site 1 being occupied 
(resp. empty) at time t. We calculate P A (N) and P B (N) as well: 

P t A (N € N) = C z t ^—^ (W\(AD + K 2 DE) N - l KD\V) = C z * A(1 - A) [A(l -p + PA)] N ~\ (32) 

1 — 2 (1 — Zj 

P t B (N € N) = P t A+B (N) - P t A (N) = C z t ~ ^ [A(l - P + /3A)]"- 1 , (33) 

P t B (0) = Q t (0)=C 2t ^. (34) 



We also set Pf(Q) = 0. Indeed Eqs. (31), (32) and (33) This result can be interpreted as follows: 



agree with the results derived i n [12] in a more compli- 
cated way. By using the result (31 1, the average number 



of particles at time t is found to be 



(Nt, 



N>1 



A+B 



C.t 



(N) 
A 



(l-z)(l-A)(l + /3A) 
"(!-") ( a < -J-) 



r(l+/3)3 



1+/3 



(a 



(35) 



(36) 



for i — > oo. 

Now we turn to the waiting time, i.e. the time that 
a new particle stays in the system. For a given number 
N of particles in the system, the probability that the 
waiting time is T € N is given by 



(A) 
(B) 



T -N 

N 
T — N 
N 



p N+1 (l 
- 1 



P) 



T-2JV 



A(T,N), (37) 



p»+\i-P) 



T-2N-1 



B(T, N). 
(38) 



Here A (resp. B) corresponds to the case where the right- 
most site is occupied (resp. empty). Note that ( x ) de- 
notes the binomial coefficient, which should not be con- 
fused with a two-dimensional column vector. The aver- 
age waiting times for given N in the cases A and B are, 
respectively, 



(T NA )= TA(T,N) = 
(Tn,b)= TB(T,N) = 



N + l 



E 

T>2AT+1 



P 

N- 



N-l, 



N. 



(39) 



(40) 



• For each particle, it takes one time step to move 
from site 2 to site 1. For N + l particles, it takes, 
in total, N time steps (resp. AT + 1 time steps) for 
the case A (resp. B). 

• For each particle, it takes i time steps in average to 
leave the system after arriving at site 1. For N + l 
particles, it takes, in total, =^±1 time steps. 

• A new particle entering the system at time t does 
not wait during time t and t + 1. Thus we have to 
subtract 1 from the above. 

Let us consider the probability W t (T) of the waiting 
time T for a particle entering the system at time t. Using 
Eqs. (pi, ((331), pi}, (137} and pH), we find 



LT/2J 



W t {T)= [A{T,N)P A {N)+B{T,N)P B (N)] 



N=0 



(41) 



where |_ - J denotes the floor function, i.e. \T/2\ = T/2 
(if T e 2N) or [T/2\ = (T - l)/2 (if T e 2N - 1). In 
the convergent phase, W t (T) converges to the stationary 
distribution of the waiting time 



lim W f (T) = lim P(l - A) (1 - P + pA) 

t—t-oo z— >1 

P — a — aft 1 1 — P + aP 



T-l 



1-P + aP 



1 



(42) 



which agrees with the result in [TU] . 

To finish this section, we investigate the average wait- 
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ing time: 



{%) = C zt ^T(l-0 + 13 Af-' 



T>1 



c z 



/3(l-z)(l-A)' 



(43) 



The order of the closest singularity z = 1 to the origin 
depends on the parameters (a, 0) [12] : 

v 1 - z 1 - a , 

» 1 Ti /3(l-z)(l-A) = /3-a-a/9 ( " < 1+^' (44) 



(l-z)t 



1 + 



lTi^(i-z)(i-A) y/3(i + /3) 

(1-z) 2 a - /? + a/3 j8 
& fli -,)(!- A) = ( a >T+^)'( 4fi ) 



and thus we have 
1 - a 



(Tt] 







!-n-n.l (a< l + /3 ) ' 



(47) 



(Tt)=2j _ a * , +o(Vt) (a=^), (48) 



7T/3(l + /3) 



1 + /3' 



, , a — /3 + a/3 , . , s , - 

(Tt) = -p -t + o(t) («> r ^g), (49) 



as i -+ oo. We note that one of the central results of 
queueing theory, Little's theorem [6], is indeed satisfied 
in the convergent phase (a < jxg) [TU] : 



a lim (T t ) = lim (N t ). 



(50) 



We also notice that, in the divergent phase and on the 
critical line (a > j+q) as well as in the convergent phase, 
there is a physically natural relation between the average 
waiting time and the average number of particles: 



J lt -(T t )~(N t ) (i->oo). 



(51) 



Here J\ t is the current of particles passing through the 
exit (outflow), which will be derived in the next section. 
Note that this relation also holds for the usual queueing 
process, see Appendix. 



IV. DENSITY AND CURRENT 

We consider the probability pj t that the site j is oc- 
cupied at time t, i.e. the density profile. The initial con- 
dition implies that pj t = for j > t and p u = oi l . The 
density profile for general j and t can be calculated as 



Pjt = Pt(lTj-i---n)+ Y Pt(lT L -i ■ ■ ■ Tj+ilTj-i ■ ■ ■ n) 

Tfc=0,l Tfc=0,l 

C z t\^—\A*{W\D(D + E)1- l \V)+ A L {W\D(D + E) L - i - 1 D(D + Ey- 1 \V) 



c z 



1 - z 
I- A 



(52) 



— -A^(W\D{D + E) s ~ l \V) + A(W\D [1 - A(T> + E)]' 1 D(D + E^-^V)} 



A^ 



(l-z)(l + /3A)- 



By definition, the particle current J lt passing through 
the exit (the right end) during t and t + 1 is given by 



Jit = 0Pi 



(53) 



The particle current Jj t through the bond between the 
sites j(> 2) and j — 1 



J 3t= J2 Pt(10^_ 2 ---Tl) 



r k =0,l 



Pt{lT L -l---T j+l WT j -2---T 1 ) 



(54) 



r k =0,l 
£>j + l 



also satisfies the relation 



J 



jt 



Pp. 



jt 



(55) 



since DE(D + E) j - 2 \V) = 0D{D + Ey-^V) . 

For the generic choice of parameters a and 0, the den- 
sity profile near the right end converges as 



Pjt 



Urn 



AJ 



i 

1+0 



fll + M (l-a) 



(l-a)/3 



0^ l+p). 



(56) 



G 



for t — > oo. Here we took the limit with the site number 
j independent of time t. In particular, for j = 1, we have 



lim p lt 



lim J u 



1 

1+/3 



g 
1+0 



( a <l|g)» 
(«<!{?). 



(57) 
(58) 



confirming the relation (51 ) 



Let us now consider rescaled density profiles in the di- 
vergent phase and on the critical line, where the average 
system length grows of order t and respectively [see 
Eqn. ( 26 )]. First we observe that the density profile pj t 
can be interpreted as the expected number of nonintcr- 
acting asymmetric random walkers at time t on site j 
since the expression ( 52 ) satisfies the equation 



Pj,t+i = apj-i,t + IPjt + Spj+i.t 



(59) 



with 7 = (1— a)(l— /3) and 5 — (1— a) (3, which is the same 
as for Q t (L) [cf. pO])]. We extend the domain j e N to 
j £ Z so that pjt can be regarded as the expected number 
of walkers with the initial condition that a random walker 
exists at each site i G Z< with probability 



1 



PiO 



1 + /3 



+ (!-<*) 



(1 - a)/3 



1-a-aP 1 

(60) 



Each walker at site j hops to its left site j + 1 with proba- 
bility a, to the right site j — I with probability (5, or stays 
at site j with probability 7 (a + 7 + 5 = I), see Fig. [3] Let 

Cj t J be the probability that the walker starting from the 
site i is in site j at time t, which is distributed around 
Vt + i as IT51 



e (i) 



f 



exp 



(j -vt- if 



2at 



(61) 



for the generic case < a < 1 and < j3 < 1. Here 
(7 = a + (5 — (a — <5) 2 , and = a — 5 which is equal to 



the velocity for the system length, see Eqn. (26). The 
density profile is expressed as 



Pjt 



Pivtjt- 

i<0 



(62) 



In the divergent phase a > (with a < 1 and < 



.. . I ~' 3 

ft < 1), noting the initial condition 
lim 



1 + 



(63) 



and the form (61), we find that the density profile with 



rescaling of the position j = xt converges as 



Pxt,t 




(0 < x < V), 
(V < x < 1). 



(64) 



3 + 



a 




f 


5 
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FIG. 3: Schematic picture of the noninteracting-random- 
walker interpretation. 
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FIG. 4: Rescaled density profile in the divergent phase. The 
parameters are chosen as (q,/3) = (4/5, 1/2). The markers x 
correspond to Eqn. (52) with t — 200, £ = j/t, and the line 



to the asymptotic form ( 64 1 



Figure [4] gives an example for the rescaled density profile 
in the divergent phase. 

On the critical line a — j^jj (0 < j3 < 1), noting the 
initial condition 



lim p i0 



1 - X 



1 + /3 



(65) 



7 



Q 




I? 0.1 
- 

Q n. 



FIG. 5: Rescaled density profile on the critical line. The 
parameters are chosen as (a,/3) = (1/3, 1/2). The markers x 
correspond to Eqn. ( 52 1 with t = 400, x — j / %/t, and the line 



to the asymptotic form ( 66 1 



FIG. 6: Rescaled density profile for (a,/?) = (4/5,1). The 
markers □ and x correspond to Eqn. ( |52[ ) with t = 60 and 
61, respectively, and x — j/t. The line corresponds to the 
asymptotic form ( 70 1. 



and the form (61) with V = 0, we find that the density 



profile ( 52 ) with the rescaling x 



converges as 




(66) 



Here erfc is the complementary error function: erfc(a;) 



= e y dy. Figure 



gives an example for the 



J2 

rescaled density profile on the critical line. 

Now we consider some special cases. When a = 1, the 
position of the leftmost particle is t by definition, and we 
have pu = I- In this case, A = z and the density profile 
becomes simply 



=C~t 



(l-z)(l + /3z) 



i - 



(67) 



In particular, the density profile observed by the leftmost 
particle is independent of time t, and exhibits oscillations. 

When /3 = 1 and a > |, another oscillation occurs. 
Since 7 = 0, the walker starting from the site i can exist 
on the site j at time t only if j — i — t G 2Z, and is 
distributed around Vt + i at time t as 



irat 



exp 



(j -Vt- if 



2at 



(68) 



We also note that the initial condition for the walker on 
site i converges as 

lim p i0 = 1 — a, lim p i0 — a. (69) 

In the limit t — > oo with the scaling j — xt, we have 

1-a (0 < x < 2a- l,j -t e 2Z + 1), 
/;,,-> ^ q (0 < x < 2a- l,j -t e 2Z), (70) 
(2a-l<x<l). 



Figure [6] gives an example for the rescaled density profile 
in this case. Also when a = | and (3 = 1, the walker 
starting from the site i can exist on the site j at time t 
only if j — i — t G 2Z, and is distributed around Vt + i 
at time t as Eqn. (68 1. Noting V = and the initial 
condition pio = 1, we find the rescaled density profile 



(66) 



For the very special case a = ft = 1, which is com- 
pletely deterministic, the site j is occupied if j < i and 
t — j is even, or empty otherwise. The density profile 
oscillates between 1 and 0. 



V. CONCLUDING REMARKS 

We have studied dynamical properties of the EQP with 
deterministic bulk hopping p = 1. We found the exact 
dynamical state in matrix product (MP) form with a two- 
dimensional representation of the matrices and vectors. 
The MP dynamical state approaches the MP stationary 



state in the convergent phase a < 



1+/9 



as t — v oo. We 



have obtained the time-dependent distributions of the 
system length L (22), the number of particles N (31) 



and the waiting time T (41), and the time-dependent 
density (52) and current (53) profiles of site j. An in- 



teresting point is that essentially they are given in the 
form 



C^izMzT (x = L,N,T,j) 



(71) 



with functions 
[cf. (24)]. 



^ and $ including the square root 



We found that the asymptotic density profile in the di- 
vergent phase (with the generic choice of parameters) is 
fiat. In contrast, the density profile for p < 1 is nontrivial 
[19]. One of the important tasks for future studies is to 
determine the form of the density profile and its depen- 
dence on the system parameters for the general case. The 
stationary state for the EQP with probabilistic hopping 



p < 1 has a matrix product form with infinite dimen- 
sional matrices [11] . This fact makes us expect that the 
MP dynamical state can be extended to the p < 1 case, 
which approaches the MP stationary state in the limit 
t — > oo as in the following diagram: 



where 



dynamical state 
of unknown form 

t — > oo I 

MP stationary state 
with ooD matrices 



MP dynamical state 
with 2D matrices 

Jt->oo 

MP stationary state 
with 2D matrices. 



In the probabilistic hopping case p < 1, however, the 
master equation cannot be simplified similar to Eqs. ([5])- 
(|7j), and the factorization ansatz Q is no longer valid 

What we have investigated here is a very basic model of 
queues with excluded- volume effect. Apart from the gen- 
eralization to p < 1 one can consider various other gen- 
eralizations of the EQP, for example multi-lane queues 
with some types of queue-changing rules. 
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Appendix A: Results on the usual queueing process 



e 



1 - (1 - a - P + 2a/3)z ■ 



(A5) 



2(1 -a)Pz 

s =Vti - (1 - a - P + 2aP)z} 2 - 4a/3(l - a) (I - P)z 2 . 

(A6) 

When a < P, the system approaches the stationary state 



lim P t (N) = lim(l - 6)6 

t— >oo z— 



A" 



P(l-a) 
P-a 



a(l - P) 
0(1 - a) 



N 



(A7) 



The average number of particles at time t is 

(N t ) = V NPAN) = C z t ^ — 

.4^ (1 - z)(l - 6) 

(a<P), 



N>1 




( "(1 


-P) 


P- 


-a 


V 


a(l— ot)t 






I (a- 


-p)t 



(A8) 



The particle current leaving the system (outflow) J t at 
time t is given by 



J t =apP t (0)+pJ2Pt(N) = C z . 



p[e + a(l-Q)] 



N>1 



1 - Z 



a (a < P), 
P (a>p), 



(A9) 



The usual queueing process is characterized by the 
number N of particles which is equal to the length of 
the system since its spatial structure is not taken into 
account. We denote the probability that the number of 
particles is N at time t by P t (N) . At each time, a parti- 
cle enters the system with probability a. When a particle 
leaves the system at time t— 1, the next particle can leave 
the system at time t with probability p. If a new par- 
ticle enters the system at time t and there is no other 
particle, this new particle can leave the system simulta- 
neously. The critical line for the usual queueing process 
is a = P, and the system is convergent or divergent if 
a < P or a > P, respectively. The dynamical solution to 
the master equation 

P i+1 (0) = [(1 - a) + ap]P t {0) + (1 - a)pP t (l), (Al) 

P t+1 (N)=(l-a)pP t (N+l) 

+ [(1 -a)(l- P) + ap]P t (N) (A2) 

+ a(l - p)P t (N - 1) (JVeN) 

with the initial condition 



P (0) = 1, P (N)=0(N€N) (A3) 



is given by 



P t (N) = C zt ^Q N , 



(A4) 



for t — > oo. Note that, in the EQP case, lim t _ i . 00 J\ t (58) 
is not equal to the exit probability P in the divergent 
phase since the rightmost site can be empty. The prob- 
ability of the waiting time T for a given number N of 
particles is (^)/3 Ar+1 (l — P) T ~ N from which we find the 
probability of the waiting time for a given time t: 



N=0 



T 



Y,P t (N)( N )P N+1 (l-P) 



\T-N 



_C«>ze> (1 _, + fl8 f. 

(1 - Z) 

The average waiting time is then 

\ z ) T>1 



(A10) 



(All) 



= C,* 



i- p + pe 
/3(i-*)(i-e) 



i-p 

P-a 



(a < P), 



a-p. 



(a>P), 

and the relation J t (T t ) ~ (N t ) (t — !• oo) holds. 
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